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Abstract. The growth of agricultural and mining activities in the Amazon has impacted 
land-use and caused significant changes in the local environmental conditions of streams. 
In the face of these changes, our study aimed at assessing how environmental changes af- 
fect Odonata larval assemblages in streams in the eastern Amazon. We hypothesized that 
habitat conditions in streams are strong predictors of Odonata larval assemblages. We 
sampled 30 headwater streams (1* through 3" order) in the eastern Amazon. We corrobo- 
rated our hypothesis that regional- and local-scale environmental changes are important 
predictors of the Odonata larval assemblage structure. These results indicate that envi- 
ronmental conditions within the stream channel are important to maintain Odonata larval 
assemblages, as they provide important resources for larval development. For new stud- 
ies, we recommend the assessment of temporal dynamics to evaluate whether these pat- 
terns are stable across time. Finally, evaluating various environmental scales of the original 
impact is extremely relevant for preventing the deterioration of or recuperating aquatic 
assemblages in Amazonian streams, considering the ongoing rapid environmental changes 
and deforestation in the region. Here we demonstrate that in-stream environmental condi- 
tions are important to assemblage structure and this must be considered in environmental 
restoration plans. 
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Introduction 


Environmental changes in streams, in limnological characteristics, and in habitat 
structure have significant effects on aquatic insect assemblages, especially on 
those that are highly sensitive to environmental changes, such as some species 
of mayflies (Ephemeroptera), odonates (Odonata), stoneflies (Plecoptera), and 
caddisflies (Trichoptera) (Castro et al., 2017; Faria et al., 2017; Ferreira et al., 
2017; Mendes et al., 2019). Any changes in the abundance or species composi- 
tion of these taxa in aquatic assemblages are of concern, because they play im- 
portant roles in nutrient cycles and energy flows (Castro et al., 2019). Odonates 
are important indicators of habitat quality due to the group’s high diversity and 
sensitivity to environmental changes (Carvalho & Nessimian, 1998; Mendes et 
al., 2015). Additionally, because they are predators, they play an important role 
in aquatic trophic webs (Corbet, 1999; Pires et al., 2020). 

Odonates have aquatic larvae and flying terrestrial adult stages. They are ex- 
tremely abundant and diverse, which allows them to occupy a wide variety of mi- 
crohabitats in lentic and lotic environments (Corbet, 1999). At their aquatic stage, 
changes in land-use (Mendes et al., 2017, 2019), habitat structure, and limno- 
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logical parameters (e.g., temperature, pH, conductivity, 
and dissolved oxygen) can result in important changes 
in assemblage properties, such as richness, species 
composition, and abundance patterns (Mendes et al., 
2017; Pires et al., 2020). These changes in the assem- 
blage structure are related to their physiological needs: 
the characteristics of thermoregulation in adults (De 
Marco et al., 2015) and the characteristics of breath- 
ing in larval instars (Corbet, 1999). Therefore, changed 
landscapes usually have a greater impact on adults (De 
Marco et al., 2015) and in-stream changes more affect 
the immature stages (Mendes et al., 2017, 2019). 

Our study aimed at assessing how environmental 
changes affect Odonata larval assemblages in streams 
in the eastern Amazon. We hypothesized that habitat 
conditions in streams are strong predictors of Odona- 
ta larval assemblages. Our hypothesis is based on the 
premise that in-stream environmental conditions are 
altered in regions that experience anthropogenic pres- 
sures and that this affects aquatic insects (Martins et 
al., 2018; Castro et al., 2019). This is due to a hierar- 
chical mechanism, where climate and dispersal affect 
metacommunities on a landscape scale, contributing to 
regional and biogeographic patterns (Brasil et al., 2019), 
whereas local habitat characteristics are important for 
the spatial organization on smaller scales (Montag et 
al., 2019; Oliveira-Junior et al., 2019). 


50.00°W 


49.00°W 48.00°W 47.00°W 






” 
ro) 
~ 
N 
9 
O 
ive) 
A 

Soil classification 2 
= EE Primary forest & 
E Pasture 
e Samples EE Exposed soil 


— Water mass Mill Secondary forest 





Habitat conditions in streams influence Odonata larval assemblages 


Material and methods 
Study area 


Our study was conducted in the Rio Capim river basin, 
in the municipality of Paragominas, in the southeast 
of the state of Para (Figure 1). The region has approxi- 
mately 150,000 ha of forested area, of which 18,000 
are pasturage, 11,000 are Permanent Protection Areas 
(PPA), 15,000 are areas where logging takes place, and 
98,000 ha are native forest. The Capim River basin has 
a drainage area of around 37,000 km?, mostly due to 
the slope of the region, and has a roughly rectangu- 
lar shape (Lima & Ponte, 2012). The forest rerpesents 
dense tropical rainforest (Almeida & Vieira, 2010) in 
a humid equatorial climate with an average rainfall of 
1,743 mm per annum. The rainy season extends from 
December through May, and the dry season from June 
through November (Alvares et al., 2013). The average 
relative humidity is around 81%, and the average an- 
nual temperature is 26.3°C (Francez et al., 2009). 
Beginning in the second half of the 20" century, the 
municipality of Paragominas has been suffering high 
rates of deforestation due to the intensification of an- 
thropogenic activities. As a result, the Rio Capim river 
basin today displays a mosaic of different land-uses, 
such as areas for agriculture and cattle farming. All 
these activities contribute to the deforestation of Bra- 
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Figure 1. Location of the 30 sampling sites in the Rio Capim river basin in the municipality of Paragominas, state of Para, Brazil. 
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zil’s Legal Amazon, which has so far been recorded as 
amounting to around 80% (Gardner et al., 2013; Olivei- 
ra-Junior et al., 2017). 


Sampling design 


The aquatic ecosystems of the Amazon region are high- 
ly variable in size and shape. Thus, we standardized 
our sampling procedure by selecting 30 streams of the 
first through third orders, according to Strahler’s clas- 
sification (1957). Surveys were carried out in the years 
of 2015 and 2017, always in the periods of less rainfall 
to reduce possible seasonal effects. In each stream, we 
defined a 150-meter stretch divided into 10 sections of 
15 meters, each separated by a cross-sectional transect 
(Figure 2) denominated as “A” through “K” from down- 
stream to upstream. 

Within each of the 10 sections, we sampled Odonata 
larvae and recorded the physical and physicochemical 
habitat parameters following the physical habitat moni- 
toring and assessment protocol of the United States 
Environmental Protection Agency (US-EPA) (Kaufmann 
et al., 1999; Peck et al., 2006). This protocol records en- 
vironmental information in different categories: avail- 
ability of shelter for aquatic communities, the structure 
of riparian vegetation, channel morphology, human im- 
pact, and water limnological variables. Evaluating en- 
vironments according to this protocol has been widely 
used in ecological studies on assemblies of aquatic in- 
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sects on and in Amazonian streams, and thus has come 
to be an important tool for assessing the effects of land 
uses in the Amazon biome (Juen et al., 2016; Mendes 
et al., 2017, 2018). The physical habitat characteriza- 
tion along each section was carried out following the 
stream assessment methodology of Peck et al. (2006) 
and Hughes and Peck (2008). This methodology has 
already been successfully applied in several studies of 
Amazonian streams (Juen et al., 2016; Mendes et al., 
2019; Montag et al., 2019; Oliveira-Junior et al., 2019; 
Dias-Silva et al., 2020). 


Biological sampling 


We furthermore subdivided the cross-sectional sec- 
tions into three five-meter segments for the biologi- 
cal sampling of individuals of which only the first two 
segments were sampled (Figure 2). Thus, we sampled 
20 segments of five meters in width in each stream. 
Sampling was carried out by sweeping the stream three 
times from the streambed towards the stream margins 
using a sieve of 18 cm of diameter and 250 um in mesh 
width (Shimano & Juen, 2016; Mendes et al., 2019). 
Sampled individuals were taken to the laboratory, pre- 
served in 85% ethanol, and identified at genus level 
using taxonomical dichotomous keys (Heckman, 2006; 
Hamada et al., 2014). To standardize the reliability of 
sample identification, only organisms that exhibited the 
development of wing pads reaching the fifth segment 


Downstream 





Figure 2. Diagram of the sampling design in each of the sampled streams in the Rio Capim river basin in the municipality of 
Paragominas, state of Pará, Brazil. We defined a 150-meter transect in each sampled stream, which was then divided into ten 
15-meter sections, each of them separated by a cross transect named “A” through “K” from downstream to upstream. 
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of the dorsal abdomen were included (Mendes et al., 
2017, 2019). Depending on the species, Odonata larvae 
may have up to 15 instars, which obstructs larval iden- 
tification at species level. Because many species do not 
have their larval instars described, identification at ge- 
nus level was a recommendation from specialists in the 
taxonomic classification of Odonata larvae that we fol- 
lowed. Although the taxonomic resolution at genus lev- 
el is currently considered the best for this group from 
the Amazon region (Mendes et al., 2017), the identifi- 
cation of some genera has remained quite challenging, 
and it is not as yet possible to distinguish, for example, 
individuals of Hetaerina from those of Mnesarete, 
which were therefore indiscriminately considered Calo- 
pterygidae. All sampled organisms were stored in the 
Zoology Collection of the Federal University of Para, 
UFPA (Belém — PA, Brazil). 


Characterization of land-uses and vegetation cover 


We used images from the Shuttle Radar Topography 
Mission (SRTM) project with a resolution of approxi- 
mately 30 meters. These images were adjusted using 
images from Google Earth (http://earth.google.com) 
and TauDem version 5.3 by means of the geoprocess- 
ing software QGis version 2.18 (QGIS 157 Development 
Team, 2017). 

The interpretation of the different types of land-uses 
and vegetation cover was performed using images from 
the Landsat 8 satellite, using the Semi-Automatic Clas- 
sification plug-in of the software QGis 2.18 (Macedo et 
al., 2014). The images were obtained from the United 
States Geological Survey (USGS) Earth Explorer project 
in 2015 and 2017. This set of images is essential, since it 
provides data on the shapes and textures of the differ- 
ent elements of the images for identifying the various 
types of vegetation cover along the streams (Castro et 
al., 2017). 

To minimize the effects of atmospheric reflectance, 
we applied a atmospheric correction process to the im- 
age set (Antunes et al., 2012). Later, the images were 
classified into primary and secondary forest areas, 
pasturage, and areas of exposed soil (where mining ac- 
tivities are concentrated). We validated the classifica- 
tion process by inspecting images from Google Earth 
(http://earth.google.com). Finally, we cropped a band 
of 30-meters wide and 600 linear meters in length of 
riparian zone from each of the stream margins, fol- 
lowing the drainage network. This band of vegetation 
(30 x 600 m) is the spatial scale widely used in aquatic 
ecology studies because it is very important for aquat- 
ic biota in Amazonian streams (Montag et al., 2019a, 
2019b; Benone et al., 2020; Leão et al., 2020). In ad- 
dition, 30 meters is the width in which Brazilian law 
requires the preservation of native vegetation on small 
streams. Therefore, our landscape matrix was com- 
posed of land use percentages such as agriculture and 
livestock, exposed soil, and native forest. These uses 
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were measured on the spatial scale explained above 
(30 x 600 m). 


Physical habitat characterization 


We measured riparian vegetation structure, limnologi- 
cal variables, channel morphology, and anthropogenic 
impact. Each section was characterized by 10 crosswise 
equidistant measuring points along the section. We 
recorded the depth of the thalweg, number of water- 
retaining objects, presence of fine sediments, slope and 
sinuosity of the channel, and type of channel. In each 
of the sections, we measured the depth and type of 
substrate through visual inspection: sand, clay, consoli- 
dated clay, leaf banks, coarse gravel, fine gravel, con- 
crete, rock formations, wood, particulate organic mat- 
ter, macrophytes, and/or algae, roots, pebbles, and silt. 

The limnological parameters temperature, pH, tur- 
bidity, dissolved oxygen, total dissolved solids (TDS), 
and electrical conductivity were measured with the 
multiparameter probe Horiba® U-50 (Ferreira et al., 
2017). Riparian vegetation structure and density were 
estimated by visually inspecting both margins within a 
band of 100 m? for the following three types of veg- 
etation: ground vegetation (< 0.5 m height), understory 
vegetation (0.5-5 m height), and canopy vegetation 
(>5 m height). 

We also assessed habitat complexity in each section 
through visual inspection of the surface covered by 
macrophytes, roots, leaves, algae, and wood fragments. 
We used a spherical densitometer to measure the per- 
centage of canopy cover over the stream main channel 
in the central portion of each cross transect. We took 
six measurements in the central portion of each tran- 
sect: left, right, center, center left, center right, center 
upstream, center downstream. We took a total of 69 
different measurements of the physical habitat in each 
sampled stream based on Kaufmann et al. (1999) (Sup- 
plementary material S4). 


Data analysis 

We performed a “forward selection” of the R software 
(R Core Team, 2020) adespatial package (Dray et al., 
2018) for the variables that best explain the variation 
in Odonata assemblages to reduce the number of habi- 
tat and landscape variables (Blanchet et al., 2008). This 
selection uses the species matrix and the matrix of en- 
vironmental predictors and was suggested by Borcard 
et al. (2018) for the selection of variables to reduce the 
residual portion in analyses such as RDA and CCA. Af- 
ter that, to avoid autocorrelation among the selected 
variables, we applied Spearman correlations to exclude 
highly correlated environmental variables (ps > 0.70). 
We kept the one with greatest importance to odonate 
assemblages for each pair of highly correlated variables. 
This prior variable selection procedure is important to 
reduce the residual variability of our final explanatory 
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model, which could be inflated by a high number of 
predictor variables (Legendre & Legendre, 2012). We 
considered each stream our sample units and the as- 
semblage data, and we summed up the abundances 
recorded in all sampled segments within the same 
stream. More details on environmental variables are 
provided in Supplementary material. 

We ran a Redundancy Analysis (RDA) to verify the 
relationship between assemblages and Odonata gen- 
era with the environmental variables selected in the 
previous step. To do so, we used these environmental 
variables as predictor and the genera composition as 
our response matrix. Before running the analysis, we 
Hellinger-transformed our abundance data (Legendre 
& Gallagher, 2001) to maximize the explanatory pow- 
er of our model (i.e., adjusted R?) (Peres-Neto et al., 
2006). We built a biplot using the RDA axes to visualize 
the association between the most important environ- 
mental predictors and the odonate assemblages. We 
also applied ANOVA to test for the linear effects of the 
environmental predictors on the main RDA axis, which 
represents the biotic matrix (Legendre & Legendre, 
2012). 

Finally, to directly test the relationship of Odonata 
communities with environmental conditions and to 
complement what was seen in the RDA, we carried out 
a Mantel test. The Mantel test was performed with a 
species composition matrix with abundance data us- 
ing the Braycurtis distance and with a matrix with en- 
vironmental variables using the Euclidean distance. All 
analyses were run in the ‘vegan’ package, ver. 2.5—-6 
(Oksanen et al., 2020), adopting a = 0.05. 


Results 
Assemblage description 


We sampled a total of 444 individuals representing 37 
genera, of which 11 (29.7%) genera and 115 (25.9%) in- 
dividuals were damselflies, and 26 (70.3%) genera and 
329 (74.1%) individuals were dragonflies. Most abun- 
dant among the nine families we sampled were Gom- 
phidae (Anisoptera), with 42.34% of the individuals, fol- 
lowed by Libellulidae (Anisoptera) and Coenagrionidae 
(Zygoptera) with 26.58% and 18.69% of the individuals 
sampled, respectively. 

The most abundant genera of the Gomphidae were 
Zonophora (Selys, 1854), with 55 individuals (29.25%), 
Progomphus (Selys, 1854), with 45 individuals (23.94%), 
Phyllogomphoides (Belle, 1970), with 37 individuals 
(19.68%), and Phyllocycla (Calvert, 1978), with 30 in- 
dividuals sampled (15.96%). The genera Gynothemis 
(Calvert in Ris, 1909) (21.19%) and Oligoclada (Karsch, 
1890) (20.34%) were the most abundant ones within 
the Libellulidae, with 25 and 24 individuals, respec- 
tively. The genera Acanthagrion (Selys, 1876) and Argia 
(Rambur, 1842) were the most abundant ones within 
the Coenagrionidae with 24 individuals (28.91%) sam- 
pled each (Supplementary material S1). 
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Table 1. “Forward selection” results for the environment pre- 
dictors related to larval odonate assemblages on streams in 
the Rio Capim river basin, in the municipality of Paragominas, 
Para, Brazil. AdjR? = adjusted R?. 








Variable Code AdjR? P 
Average intermediate woody XCMW 0.072 <0.001 
canopy 
Temperature TEMP 0.117 0.003 
Average natural shelter (leaf banks XFC LIF 0.150 0.013 
and live roots) = 
Volume of wood in the streambed V3W_150 0.180 0.020 
Percentage of algae PCT_AL 0.208 0.022 
Proportion of shelter at excavated PFC_UCB 0.208 0.050 
margins 
P ti fh i t 

ropor on of human impact per ac 0260 0.011 
section E 





Association between assemblages and the environ- 
mental predictors 


The “forward selection” procedure selected seven local 
environmental variables that apparently were related 
to the odonate assemblages: average intermediate 
woody canopy, average natural shelter, percentage of 
algae, proportion of human impact per section, propor- 
tion of shelter at excavated margin, temperature, and 
volume of wood in the streambed (Table 1). 

When we looked into the association between en- 
vironmental predictors and odonate assemblages, we 
found that they explained 25% of the variation con- 
tained in the first two RDA axes (ANOVA for test linearity 
relationship between environmental and biotic matrix; 
F = 2.45; df = 7; p < 0.001). Both axes were negatively 
affected by temperature and the volume of wood in the 
streambed and positively affected by natural shelters 
(leaf banks and live roots) (Supplementary material S2). 
According to the Mantel test, used as a complement 
to the RDA, the environmental variables together ex- 
plained 24% (Ruane: = 0-245; p = 0.001) of the distribu- 
tion of Odonata communities in the streams analyzed. 

The RDA graphical representation of the association 
between the Odonata larvae assemblage and local and 
regional environmental variables showed, in the first 
axis, a gradient from streams with larger amounts of 
woody vegetation (negatively related) to those with 
higher variation in algae percentage (positively related) 
(Figure 3). The genera Aeschnosoma and Phyllocycla 
were closely related to streams with a higher variabil- 
ity in the amount of algae. However, the second axis 
showed a gradient from streams with higher water 
temperatures (negatively related) to streams with a 
higher variability in natural shelters, such as leaf banks 
and roots (positively related). The genus Zonophora 
was more closely related to streams with a larger num- 
ber of natural shelters and a higher percentage of algae 
(Figure 3). 
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Figure 3. Redundancy analysis (RDA) plot showing the associations between local habitat predictors and landscape predictors 
with the Odonata larval assemblage. Different colors represent the different land-use types: Dark green: forest fragments; 
Light green: secondary forest; Yellow: pasturage; Red: exposed soil. 


Discussion 


The variation in Odonata larval assemblages is explained 
by the composition of the local habitat, supporting our 
hypothesis. This is a recurring pattern in studies that in- 
vestigate the effects of environmental change on fish 
(Montag et al., 2019) and aquatic insect assemblages in 
tropical streams, such as mayfly, stonefly, and caddisfly 
assemblages (Siegloch et al., 2017), and adult odonate 
assemblages (Oliveira-Junior et al., 2019). Although 
changes in the landscape affect in-stream conditions 
(Leal et al., 2018), all the evidence cited above, togeth- 
er with our results, suggest that changes at small spatial 
scales (in-stream) are the ones that most directly affect 
aquatic diversity. 

All environmental predictors that were important in 
explaining assemblage variability were related to the 
stream marginal vegetation cover (canopy cover, hanging 
vegetation, leaf banks and roots, and wood fragments 
within the channel). Other predictors may be indirectly 
affected by deforestation, such as water temperature, 
which can increase due to intensified insolation (Carval- 
ho et al., 2013), decrease in thalweg depth due to higher 
sediment deposition (Wood & Armitage, 1997; Couceiro 
et al., 2010), increase in stream margin excavation due 
to the lack of protecting vegetation (Oliveira-Junior et 
al., 2017), and human impact, which is most severe in 
the form of deforestation (Sonter et al., 2017). 
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The local environmental changes caused by the al- 
teration of landscapes due to multiple land-uses were 
important predictors of the variation in Odonata larval 
assemblages (Juen et al., 2014; Garcia-Garcia et al., 
2017; Mendes et al., 2018). A possible explanation for 
this finding is that adult odonates are dependent on the 
physical habitat structure, such as canopy cover and 
hanging vegetation (Vianna & De Marco, 2012; Mendes 
et al., 2017). Studies that considered both larval and 
adult odonate assemblages in Brazilian streams found 
congruences in the spatial distribution of these two life 
stages, which varied from > 50 to > 90% (Valente-Neto 
et al., 2016; Mendes et al., 2017). Therefore, it is ex- 
pected that environmental predictors related to light 
exposure, such as canopy cover (De Marco et al., 2015), 
and the presence of perches, such as hanging vegeta- 
tion (Rodrigues et al., 2019), which mostly affect adults, 
would also indirectly affect the spatial distribution of 
their larvae. 

All taxa with the strongest relation to environmental 
variables were dragonflies (see Figure 3 and Supple- 
mentary material S3). This group is more tolerant of en- 
vironmental variation and potentially more inclined to 
colonizing habitats that are distant from those where 
they emerged (Corbet, 1999). Thus, they will benefit 
from the presence of open areas in that it facilitates 
higher foraging activity and offers a higher potential 
for dispersal and the colonization of new habitats for 
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larval development (De Marco & Latini, 1998; Corbet, 
1999). The removal of marginal vegetation also caused 
the transport of allochthonous debris and the growth of 
macrophytes and algae, which can favor the occurrence 
of certain dragonfly genera (Mendes et al., 2019). For 
example, the Zonophora and Progomphus individuals 
were positively affected by the extent of human impact 
on marginal vegetation. We believe this happened be- 
cause these organisms are fossorial and have a predi- 
lection for substrates with higher sediment deposition 
(Corbet, 1999; Carvalho & Nessimian, 1998). Therefore, 
because they have a narrower tolerance response to en- 
vironmental changes, some dragonflies exhibit a lower 
dependence on habitat integrity for both colonization 
and development on both regional and local scales. 
Our results demonstrate that local-scale habitat fil- 
ters can affect the structure of Odonata larval assem- 
blages along the microbasin. We believe that the altera- 
tion of environments may be a determinant for the con- 
trasting results we obtained regarding the numbers of 
genera and individuals of the two odonate suborders. 
Finally, we suggest more studies that consider different 
spatial scales (local and regional) be conducted in re- 
gions that are subject to multiple land-use types. Also, 
if possible, we suggest measuring the temporal dynam- 
ics of those landscapes and the behavior and prefer- 
ences of Odonata larvae, as they are of fundamental 
importance for filling the knowledge gaps concerning 
the effects of landscape changes on the environmental 
conditions of streams and their biological diversity. 
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Supplementary material S1. List of families and genera with the val- 
ues of abundance (N) and relative frequency (RF) of conformity 
of the suborders Anisoptera and Zygoptera recorded in the 30 
streams sampled in the Capim river basin, in the municipality of 
Paragominas, Pará, Brazil. 


Supplementary material S2. Loadings of the Redundancy Analysis 
with the local and regional predictors for the Odonata larvae as- 
semblage structure in streams sampled in the Rio Capim river 
basin, on the municipality of Paragominas, Pará, Brazil. 


Supplementary material S3. Results of RDA correlations for compo- 
sition of Odonata larvae of streams sampled in the Rio Capim 
river basin, Paragominas, Pará, Brazil. Bolds values are correla- 
tions higher than 0.1 and -0.1. 


Supplementary material S4. Local habitat structure variables ob- 
tained from the protocol for the 30 streams sampled in the Rio 
Capim river basin, municipality of Paragominas, Pará, Brazil. 
Code, Mean and Standard deviation are informed. 
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